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Abstract 

We derived explicit symbolic expressions for the first, second, and 
third Lyapunov coefficients of the complex focus of a planar system 
modelling activity of a neural network. The analysis of these expres- 
sions allowed us to obtain new results about the number and location 
of limit cycles in the model. 

Introduction. In this paper we study the mathematical model which de- 
scribe processes in a neural network consisting from two neurons. The model 
was suggested in [1] and is defined as a planar differential system 

ill = -Ml + qn(p{ui) - qi2U2 + ei, 
U2 = -U2 + q2iV{ui) +62, 

where gn, gi2, q2i > 0, ei, 62 G R are parameters, and 



1 + e-*"i ■ 

In [1] it was fulfilled a study of model (1) with the help of the bifurcation 
theory and it was listed possible types of dynamical behavior of the neu- 
ral network; with analytic and numerical techniques it was constructed the 
bifurcation diagram of system (1) for some parameter values. 

In [2] it was proved by analytic methods the correctness of the bifurcation 
diagram for some specific cases. 

*The work is supported by the Russian Foundation for Fundamental Research (Grant 
No. 99-01-00574) 
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In [2] it was demonstrated that the change of variables 



U2 — 62 
u — Ui,V — , 

forms system (1) to the system 

ii — —u + a(fi{u) — bv + c, 
i) — —V + (p{u), 

which involves three parameters only 

a = gii, & = gi2g2i, c = ei - ^1262, a, 6 > 0, c e R. 

We intend to study system (2). 

Using our results from [3], we derived analytic formulas for the first, 
second, and third Lyapunov coefficients of the complex focus of system (2). 
A sophisticated treatment of the formulas allows us to assert that for some 
parameter values system (2) has three concentric limit cycles around a unique 
steady state, two of these cycles being stable. Such a phase portrait is absent 
in the list of possible phase portraits of (2) from [1,2] although on our mind 
it is interesting from a physical point of view. 

1. Let us to note the following property of system (2). If u{t),v{t) is the 
solution of the system then ui{t) — —u{t),vi{t) — 1 — v{t) is the solution of 
the system 

ii = —u + aipiu) — bv — c + b — a, 
V = —V + (p{u), 

This fact means that if we have a phase portrait of (2) for parameter values 
a — ai,b — bi,c — ci then there exists a symmetric with respect to the 
point u — 0,v — 1/2 phase portrait of the system for parameter values 
a = ai,b = bi,c = C2 = — ai+61— Ci. In the parameter space points (ai, bi, Ci) 
and (ai,&i,C2) lie on different sides from the hyperplane a — b + 2 = (on 
equal distances). So, it is sufficiently to study system (2) for parameter values 
from one of a half-space, for example, from the half-space a — b + 2c> 0. 

The mentioned property is obvious from the bifurcation diagrams con- 
structed with numerical methods in [2]. 

The authors of [1,2] solved all problems of the number, location, and types 
of steady states of system (2). The system may have between one and three 
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steady states in the finite part of a pliase plane. If the system has three 
steady states then one of them is always a saddle. Two remaining steady 
states are nodes or foci (stable or unstable). We distinguish between them a 
left and right steady state. 

The authors of [4] introduced the following method to describe the phase 
portraits of system (2). The symbol s (u) denotes a stable (unstable) equi- 
librium; the symbol 5* (f/), a stable (unstable) limit cycle. The subscript 1 
refers to the symbol corresponding to the left equilibrium or to a limit cycle 
around it; the subscript 2, to the symbol corresponding to the right equi- 
librium or to a limit cycle around it. A symbol has no subscripts when the 
steady state is unique. The symbols S and U without subscripts refer also 
to limit cycles around all three steady states. If there are several cycles with 
the same subscript or without subscript, then they are listed in order from 
inside. 

The authors of [1,2] detected following phase portraits of system (2).[| 
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As a result of our study, we detect one more phase portrait of system (2) 
which is described by the sequence 

15. uSUS (4) 

and corresponds to the situation when in the system two stable limit cycles 
coexist which are separated with an unstable cycle. 

The detection is due to a more comprehensive analysis of the Lyapunov 
coefficients of the complex focus of system (2). The next section is devoted 
to this analysis. 

^The subscript 1 for the number of a phase portrait imphes that a corresponding 
portrait is symmetric to the portrait with the same number without the subscript. 
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2. In [3] we derived symbolic expressions for the first, second, and third 
Lyapunov coefficients of the complex focus of the system 



in terms of the coefficients aij. 

The expressions are rather unwieldy and in general case are not amenable 
to theoretical study but in some special cases they may be simplified notice- 
ably. In particular, if we have a Lienard system with the complex focus in 
the origin 

^ = D-, /g\ 
y = p{x) + yq{x) = pix+p2X^ + . . . + y(qix + q2x'^ + ■■■), 

then the expressions for the first three Lyapunov coefficients with accuracy 
to a positive factor coincide with ones 

k = 5(^2^3 - P3?2) + 3(p4gi - Pm), /gx 
14p2(P2g4 -P4g2) + 21pi(p3?4 -P4g3) + 35pi(p5g2 -P2g5) + 

+I5pi{piq6-p6qi). 
Let (uq, vq) be a steady state of system (2). Then we have 



x = y, 

y = aiox + aoiy + a2ox'^ + auxy + ao2|/^ + ■ ■ ■ 



Mo - (a - b)ip{uo) -c = 0,vo = ^p{uo). 



Let us make the following change of variables in system (2) 
X — u — uq, y — —u + a(fi{u) — bv + c. 



After it the steady state {uo,Vo) transfers to the origin. 

We obtain a Lienard system which is equivalent to system (2) 



x = y, 

y ^p{x)+yq{x), 



(7) 



where 



p{x) = c + (a — 6)/(l + exp(— 4mo — 4,x)) — uo — x, 

q{x) = -2 + (4aexp(-4Mo - 4a;)/(l + exp(-4Mo - 4a;))^. 



(8) 
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The origin will be the complex focus of system (7), (8) if p(0) = q{0) = 
0,p'(0) < 0, that is under the fulfillment of the conditions 

c + (a - b)/{l + exp(-4Mo)) -uo = 0, 

-2 + (4aexp(-4Mo)/(l + exp(-4Mo))^ = 0, 

-1 + 4(a - b) exp(4Mo))/(l + exp(4Mo))^ < 0, 

which may be written as 

c = (6 - a) /{I + exp(-4Mo)) + uq, 

a = (1 + exp(4iio))V(2exp(4xio)), (9) 
-1 + 4(a - b) exp(4iio))/(l + exp(4iio))' < 0. 

Under conditions (9) system (7), (8) takes form (5). 

Evaluating needed coefficients pi, qi and substituting them to (6), we ob- 
tain the expressions for the Lyapunov coefficients of system (7), (8). 

After some transformations we obtain that the Lyapunov coefficients with 
accuracy to a positive factor are equal to the expressions 

1-2 exp(4Mo) - (6 - 8a + 86) exp(8Mo) - 2 exp(12Mo) + 
+ exp(16Mo), 

^2 = 3-72 exp(4Mo) + (45 + 8 (a - 6)) exp(8Mo) + 

(240 - 368(a - b)) exp(12Mo) + (45 + 8(a - b)) exp(16iio)+ 
-72exp(20Mo) + 3exp(24Mo), 
= 1 - 4(29 + a - 6) exp(4Mo) + (717 + 160(a - b)) exp(8Mo) + 

+ 16(102 - 205(a - 6) - 6(a - bf) exp(12Mo)- (10) 

-2(903 - 1904(a - 6) - 384(a - bf) exp(16Mo) + 

+8(651 + 1813(a - b) - 1608(a - bf) exp(20Mo)- 

-2(903 - 1904(a - 6) - 364(a - 6)^) cxp(24?io) + 

+ 16(102 - 205(a - b) - 6(a - bf) exp(28Mo) + 

+ (717 + 160(a - b)) exp(32Mo) - 4(29 + a - b) exp(36Mo) + 

+ exp(40Mo)- 

Conditions (9) define a two-dimensoinal bifurcation manifold in the pa- 
rameter space. This manifold corresponds to the codim 1 Andronov-Hopf 
bifurcation. The condition li — from (10) defines in this manifolds a one- 
dimensional bifurcation codim 2 manifold corresponding to the degenerate 
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Andronov-Hopf bifurcation (the Bautin bifurcation, the Takens bifurcation, 
the vanishing of the first Lyapunov coefficient). 

The last manifold is piecewise smooth and it may be defined as follows 

a = (1 + ^)2/(2^9), 

c= (l-3i?-3T?2 + i?3 + 2i?lni?)/(8t?), 
< < +00, t?^ 1. 

Coordinates {uq, vq) of a corresponding unique equilibrium may be ob- 
tained from equations ex.p{Auo) = 'd,Vo = (p{uo). 

If parameters a,b,c are evaluated accordingly to (11) then system (7), 
(8) (and then system (2) also) has the complex focus of the multiplicity 2. 
Stability of the focus is defined by the sign of the expression I2 which is the 
second Lyapunov coefficient I2 from (10) evaluating accordingly to the (11) 
and to the equality -3 — exp{Auo). 

k = 2(1 + - IM + 6^^ - IM^ + ^^). 

The equation ^2 = has two positive roots ■j?!, ■j?2 

1^1 = {7 + 3^5 + v^7 l5 + v^ )/2 ^ 13.6349, 
^2 = {7 + 3^5 - \/6^ 15 + V5) /2 ft; .0733414, 

to which correspond the following parameter values oi ^ 7.8541, 61 26.9164, 
ci ^ 18.4129 and 02 = ai, 62 = C2 ~ 0.64937, ai-6i + 2i > 0, 02-62 + 22 < 
0. 

The third Lyapunov coefficient is negativ for (ai, bi, Ci) and (02, &2, C2). 

For such parameters system (2) has a unique steady state which is the 
complex focus of multiplicity 3. 

It is well known that in the parameter space to the point corresponding 
to the vanishing of the second Lyapunov coefficient is adjacent a region of pa- 
rameter values for which the system has three concentric limit cycles around 
the equilibrium. 

Thus, we proved that there are parameter values a, b, c, for which system 

(2) has three limit cycles surrouding a unique steady state (an unstable fo- 
cus) . Two from these cycles are stable and they are separated by an unstable 
cycle. 
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3. Fig.l gives a part of a complete bifurcation diagram of system (2) for 
a fixed value of parameter a > Ui. The rest of the diagram may be easily 
restored from the about mentioned reasons of the symmetry. The diagram is 
fairly typical. In particular, it has many common features with the diagram 
given in [1,2] for a < 04. We use the notations introduced in [1,2] and restrict 
ourselves to only brief remarks. More detailed description of the mentioned 
bifurcations may be found in [6,7], for example. 

The curve srii corresponds to the availability in the system of a a double 
steady state which by an appropriate parameter perturbation splits into a 
saddle and a node which lies lefter then a saddle. On the curve sni there 
is a point tbi, which corresponds to the Bogdanov-Takens bifurcation and 
divides the curve into two segments with stable and unstable saddle-nodes. 

The line hi {h2) corresponds to the Andronov-Hopf bifurcation of the 
left (right) steady state if parameters are inside the region of multiplicity 
of steady states or of a unique steady state. As we noted above these lines 
are defined by formulas (9). The point dhi on the line hi corresponds to 
the degenerate Andronov-Hopf bifurcation and divides the line into two seg- 
ments corresponding to a soft or hard loss in stability. The second Lyapunov 
coefficient is positive at the point dhi. Let us note that this coefficient is 
negative at the point dhi on the diagram presented in [1,2] for a < ai. 

The curves sli and s/2 correspond to separatrix loops around left and 
right steady states. The curve 5/4 corresponds to a separatrix loop around 
both steady states. 

The curve snpo corresponds to a double cycle which is stable from out- 
side. The curve snpoi corresponds to a double cycle which is stable from 
inside. The curves snpo and snpoi coalesce and form a cusp at the point 
corresponding to a triple cycle of system (2). 

The listed curves divides the parameter space into regions within which 
the system has distinct rough phase portraits describing in (3), (4). 

As usual, we give only the scalclcss relative positions of the bifurcation 
sets. The drawing of these sets to real scale involves technical dificulties. 

The local bifurcation sets srii, hi, tbi, dhi can be described by explicit for- 
mulas (see (9), (11), for example). Therefore, construction of these sets 
involves no difficulties. As to sets of nonlocal bifurcations sLi, snpo, snpoi, 
in constructing them we used the results of the bifurcation theory, and the 
results of the qualitative theory of ordinary differential equations, and the 
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results of simulations. 

A considerable difference of our bifurcation diagram from the diagram in 
[1,2] (for a < ai) is the presence of region 15 and curve snpoi whose existence 
we proved above. 

In Fig. 2 is shown the phase portrait of system (2) for parameter values 
a = 16., b = 130., c = 111.165, corresponding to region 15.0 




0.5 1 1.5 2 2.5 3 3.5 



Fig.2 

Hence our analysis of formulas for the Lyapunov coefficients for the system 
describing the dynamics of a neural network consisting from two sells enables 

^To confirm the existence of three hmit cycles in system (2) for the specified parameter 
values we studied the proper Poincare map numerically and established that it has three 
fixed points. 
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us to supplement results of other authors. In particular, to earlier descrided 
phase portraits we add one more portrait with three limit cycles around an 
unstable steady state. Studing a complete bifurcation diagram, we specify 
the position of this portrait in the parameter space. 

The results of [1,2] and our results give grounds to advocate that all 
possible phase portraits of system (2) are included in the join of lists (3), (4). 
4. As we noted above, the bifurcation diagram from [1,2] is fairly typical and 
occurs during a study of many mathematical models, see [8,9] for references. 
In particular, a simular diagram was consnructed by us in [10] where we 
investigated a planar system modelling a nerve conduction in the squid giant 
axon. 

Acknowledgmens. We thank Professor Fotios Giannakopoulos for ma- 
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